MATH 18.152 COURSE NOTES - CLASS MEETING # 9 



18.152 Introduction to PDEs, Fall 2011 Professor: Jared Speck 

Class Meeting # 9: Poisson's Formula, Harnack's Inequality, and Liouville's 

Theorem 



I. Representation Formula for Solutions to Poisson's Equation 
We now derive our main representation formula for solution's to Poisson's equation on a domain 

tt. 

Theorem 1.1 (Representation formula for solutions to the boundary value Poisson 
equation). Let Q be a domain with a smooth boundary, and assume that f e C 2 (Q) and g e C(dQ). 
Then the unique solution u e C 2 (Q) n C(Q) to 

(1.0.1) Au(x) = f(x), xeflcl", 

u(x) - g{x), x e dQ. 

can be represented as 

(1.0.2) u{x) = f f(y)G(x,y)d n y + f g(a) V ^ {a) G(x,a) da, 

V 

Poisson kernel 

where G(x,y) is the Green function for Q. 

Proof. Applying the Representation formula for u Proposition, we have that 

(1.0.3) u{x)= f §{x-y)f{y)d n y- f <$>(x - a)V M {a) u{a) da + f g{a)S7 ff {a Mx - a) da. 
Recall also that 

(1.0.4) G(x,y) = $(x-y)-<f)(x,y), 

where 

(1.0.5) A y( j)(x,y) = 0, xeQ, 

and 

(1.0.6) G(x,a) = when x e £7 and a e dVt. 

The expression (1.0.3) is not very useful since don't know the value of V N( a )U(a) along dQ. To 



fix this, we will use Green's identity. Applying Green's identity to the functions u(y) and <j)(x,y), 
and recalling that A y (f)(x,y) = for each fixed x e Q, we have that 
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Au(y) 



(1.0.7) 0= f (f)(x,y) f(y) d n y - [ 0(x,cx)V #u(<t) da + [ g(a) V^0(x, a) da. 

Subtracting ( 1.0.7[ ) from ( |1.0.3 ), and using (1.0.4), we deduce the formula ( 1.0.2[ ). 



□ 



2. Poisson's Formula 



Let's compute the Green function G(x,y) and Poisson kernel P(x,a) d = S/^G(x, a) from (1.0.2) 



in the case that Q = B R (0) c is a ball of radius R centered at the origin. We'll use a technique 
called the method of images that works for special domains. 

Warning 2.0.1. Brace yourself for a bunch of tedious computations that at the end of the day will 
lead to a very nice expression. 

The basic idea is to hope that <p(x,y) from the decomposition G(x,y) = &(x-y) -<p(x,y), where 
4>(x,y) is viewed as a function of x that depends on the parameter y, is equal to the Newtonian 
potential generated by some "imaginary charge" q placed at a point x* e B^(0). To ensure that 
G(x, a) = when a e dBft(0), q and x* have to be chosen so that along the boundary {y e M. 3 \ \y\ = 
R}, <p(x,y) = 47r ^_ y | - In a nutshell, we guess that 



(2.0.8) 



G(x,y) 



A7r\x - y\ 4ir\x* - y\ 

V- J 

V 

<f>(x,y)? 



and we try to solve for q and x* so that G(x,y) vanishes when \y\ = R. 

Remark 2.0.1. Note that A y g^— = 0, which is one of the conditions necessary for constructing 
G(x,y). 

By the definition of G(x,y), we must have G(x,y) = when \y\ = R, which implies that 



(2.0.9) 

Simple algebra then leads to 



4:7r\x - y\ 4tt\x* - y\ 



(2.0.10) 



\x* -y\ 2 = q 2 \x-yf 



When \y\ = R, we use (2.0.10) to compute that 



(2.0.11) 



*|2 



2x*-y + R 2 = \x* - y\ 2 = q 2 \x - y\ 2 = q 2 (\x\ 2 -2x-y + R 2 ), 



where • denotes the Euclidean dot product. Then performing simple algebra, it follows from (2.0.11 ) 
that 



(2.0.12) 



\x 



*|2 



+ R 2 - q 2 {R 2 + \x\ 2 ) = 2y ■ (x* - q 2 x). 



Now since the left-hand side of (2.0.12) does not depend on y, it must be the case that the 
right-hand side is always 0. This implies that x* = q 2 x, and also leads to the equation 
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(2.0.13) 



q A \x\ 2 -q 2 (R 2 + \x\ 2 ) + R 2 = 0. 



Solving (2.0.13) for q, we finally have that 



(2.0.14) 

(2.0.15) 
Therefore, 

(2.0.16) 
(2.0.17) 



R 
\x\ 

X = r—rxX. 

\x\ 2 



<f)(x,y) 



my) 



R 



47T \ X \\&L X -y\ 



1 



AnR 



Next, using (2.0.8), we have 
(2.0.18) 



where we took a limit as x -> in (2.0.16) to derive (2.0.17) 

R 



G(x,y) = --^ r+ ' 



An\x-y\ 4tt \xl\0zx - y\ 



(2.0.19) 



G(0,y) 



1 1 
+ 



An\y\ A-kR 
For future use, we also compute that 



(2.0.20) 



V y G(x,y) 



x-y 1 R x* -y 



Air\x-y\ 3 Air \x\ \x* - y\ z 



Now when a e dB R (0), (2.0.10) and (2.0.14) imply that 



(2.0.21) 



R 



\x - a\ 



\x - cr 



x 



Therefore, using (2.0.20) and (2.0.21), we compute that 



(2.0.22) 



V a G{x 1 a) 



x - a 



1 \x\ 2 x* - a 



xtO, 



Aty\x - cr\ 3 Air R? \x - a\ 3 

z h-K) 

4vr|x-cr| 3 V R? ) 



R 

x - a 1 \x\ 2 ~ ° 
+ 



An\x - cr| 3 Air R? \x - a\ 



Using (2.0.22) and the fact that N(cr) = icr, we deduce 



(2.0.23) 



AtyR \x - cr| 3 
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Remark 2.0.2. If the ball were centered at the point p e IR 3 instead of the origin, then the formula 



(2.0.23) would be replaced with 



(2.0.24) 



N(a) 



R 2 - \x - p\ 2 1 
AttR \x - o~\ 



Let's summarize this by stating a lemma. 
Lemma 2.0.1. The Green function for a ball B R {p) c 



is 



(2.0.25a) 



(2.0.25b) 



G(x,y) = -—^ r + J 



R 



An\x-y\ ^\x-p\\^{x-p)-{y-p) 



G(p,y) = -^ : r+ ' 



A7i\y-p\ AirR 
Furthermore, if x e B R (p) and a e dBji(p), then 



x +p^ 



(2.0.25c) 



Vjv ((T) G(z,0-) 



R? - \x - p\ 2 1 



AnR 



\x - a 



3 ' 



We can now easily derive a representation formula for solutions to the Laplace equation on a ball. 

Theorem 2.1 (Poisson's formula). Let B R {p) c IR 3 be a ball of radius R centered at p = 
(p 1 ,^ 2 ,^ 3 ), and let x = (x 1 , x 2 , x 3 ) denote a point in IR 3 . Let g € C(dB R (p)). Then the unique 
solution u € C 2 (B R (p)) n C{B R {p)) of the PDE 



(2.0.26) 

can be represented using the Poisson formula: 



Au(x) = 0, x e B R (p), 
u(x)=g(x), xtdB R (p), 



(2.0.27) 



u(x) 



R 2 - \x -pp 



L 



9 Br (p) |x - cr| : 



dcr. 



Remark 2.0.3. In n dimensions, the formula (2.0.27) gets replaced with 



(2.0.28) 



u(x) 



R 2 - \x -p\ 2 

UJr,R 



9{°) 



dB R (p) \x - a\ n 



da, 



where as usual, u n is the surface area of the unit ball in IR n . 

Proof. The identity (2.0.27) follows immediately from Theorem 1 1 . 1 1 and Lemma 2.0.1 



□ 



3. Harnack's inequality 

We will now use some of our tools to prove a famous inequality for Harmonic functions. The 
theorem provides some estimates that place limitations on how slow/fast harmonic functions are 
allowed to grow. 
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Theorem 3.1 (Harnack's inequality). Let 5r(0) c M. n be the ball of radius R centered at the 



origin, and let u e C 2 (B R (0)) nC(B R (0)) be the unique solution to (2.0.26). Assume that u is 
non-negative on B R (0). Then for any x e Br(0), we have that 



R n - 2 (R-\x\) 



R n - 2 (R+\x\) 



Proof. We'll do the proof for n = 3. The basic idea is to combine the Poisson representation formula 



with simple inequalities and the mean value property. By Theorem 2.1, we have that 

R 2 - \x\ 2 



(3.0.30) 



u(x) 



Jdi 



da. 



AirR JdB R (o) \x - a\ 

By the triangle inequality, for a e dB R (0) (i.e. \a\ = R), we have that \x\ - R < \x - a\ < \x\ + R. 



Applying the first inequality to (|3.0.30|), and using the non-negativity of g, we deduce that 

g(a)da. 



(3.0.31) 



u(x) < 



R + 


X 


1 


R 2 - 


X 


2 AttR 



Now recall that by the mean value property, we have that 



(3.0.32) 



^(0) = a ~fvj f g{a)da. 

AnR 2 JdB R (0) 



Thus, combining (3.0.31) and (3.0.32|), we have that 
(3.0.33) 



, . R(R + \x\) , n . 



(R-\x\y 



which implies one of the inequalities in (3.0.29). The other one can be proved similarly using the 
remaining triangle inequality. 

□ 

We now prove a famous consequence of Harnack's inequality. The statement is also often proved 
in introductory courses in complex analysis, and it plays a central role in some proofs of the 
fundamental theorem of algebra. 

Corollary 3.0.2 (Liouville's theorem). Suppose that u e C 2 (IR n ) is harmonic on IR n . Assume 
that there exists a constant M such that u(x) > M for all x £ M. n , or such that u(x) < M for all 
x € ]R n . Then u is a constant-valued function. 

Proof. We first consider the case that u(x) > M. Let v d = u + \M\. Observe that v > is harmonic 
and verifies the hypotheses of Theorem 3.1[ Thus, by (3.0.29), if x e M n and R is sufficiently large, 
we have that 



(3.0.34) 



R n - 2 (R-\x\) /n . , . R n - 2 (R+\x\) , n . 
/ ^ ■ is- 1 i ^(0) ^ v ( x ) ^ ttz , ; ' v(0). 



Allowing R -> oo in ( |3.0.34 ), we conclude that v(x) = v(0). Thus, v is a constant-valued function 
(and therefore u is too). 
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To handle the case u(x) < M, we simply consider the function w(x) d = -u(x) + \M\ in place of 
v(x), and we argue as above. 

□ 
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